Gap generation in the XXZ model in a transverse magnetic field 
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The ground state phase diagram of the ID XXZ model in a transverse magnetic field is obtained. 
It consists of the gapped phases with different types of long range order (LRO) and critical hnes at 
which the gap and the LRO vanish. Using scaling estimations and a mean-field approach as well as 
numerical results we found critical indices of the gap and the LRO in the vicinity of critical lines. 



75.10.Jm - Quantized spin models 

The study of the ID spin-1/2 XXZ model in a trans- 
verse magnetic field has been drawn much attention last 
years. The Hamiltonian of this model is: 



qy qy 



(1) 



The spectrum of the XXZ model for — 1 < A < 1 is 
gapless. When the transverse magnetic field is applied 
a gap in the excitation spectrum seems to open up. It 
is supposed [1] that this effect can explain the peculiar- 
ity of low temperature specific heat in Yh^As^ [2]. The 
magnetic properties of this compound is described by the 
XXZ Hamiltonian with A « 0.98 and it was shown [I] 
that the magnetic field in an easy plain induces a gap in 
the spectrum leading to a dramatic decrease of the linear 
term in the specific heat. 

At ft, = the model (I) is the well-known XXZ model. 
In the Ising-like region A > 1 the ground state of the 
XXZ model has the Neel long-range order (LRO) along 
the Z axis and there is a gap in excitation spectrum. 
In the region — 1 < A < 1 system is in the so-called 
spin-liquid phase with a power-low decay of correlations. 
Finally, for A < —1 the ground state is the classical 
ferromagnet with the gap above the ferromagnetic state. 

At ft ^ the total is not conserving and the model 
(1) is not integrable, except some special cases: A = 1 
and A — > ±cx). In addition, there is a 'classical' line 
ftci(A) — \/2{l -I- A), where the quantum fluctuations of 
XXZ model are compensated by the transverse field and 
the exact ground state of (1) is a classical one [3]. The 
excited states on the classical line are generally unknown 
(except some of them [4]), though it is assumed that the 
spectrum is gapped. 

In the limits A ±oo the model (1) reduces to the 
ID Ising model in the transverse field (ITF), for which 
the phase transition occurs at he — |A|/2. At this field 
the gap is closed and the LRO in Z direction vanishes. 

It was shown [5] that the phase transition of this type 
takes place for any A > 0. One can expect also that such 
a transition exists for any finite A at some critical value 
ft — ftc(A) and there is the transition line connecting 
two limiting points A — > ±oo. Besides, there are other 
transition lines characterizing by vanishing of both the 
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FIG. 1. Phase diagram of the model (1). The thick solid 
lines denote the critical lines, thin solid line is the 'classical' 
line, and dashed line denotes the line fei(A) (see below). 



gap and the LRO. These lines are: ft = 0, |A| < 1; 
A = 1, ft < 2; A = -1, ft < ftc(-l)- However, the 
critical properties in the vicinity of these transition lines 
are not known yet. 

Thus, we expect that the phase diagram of the model 
(1) (on (A, ft) plane) has a form shown on Fig.l. It 
contains four regions corresponding to different phases 
and separated by the transition lines at which the gap 
vanishes. Each phase is characterized by its own type of 
the LRO: the Neel order along the Z axis in the region 

(1) ; the ferromagnetic order along the Z axis in the region 

(2) ; the Neel order along the Y axis in the region (3); and 
in the region (4) there is no LRO except magnetization 
along the field direction X (which, certainly, exists in all 
above regions). 

In this paper we investigate the behavior of the gap 
and the LRO near the transition (critical) lines. We are 
interested in the critical exponents along these lines. 

The line h = 0, |A| < 1. 

Low-energy properties of the XXZ model are de- 
scribed by a free boson field theory. Therefore, to 
study the behavior of the system near the line ft = 0, 
I A| < 1, we use conformal estimations of small perturba- 
tion ft ^ 1. 

The time-dependent correlation functions of the XXZ 
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chain show the power- law decay at |A| < 1 and have the 
asymptotic form [6] 



{S^{x,t)S^ {0,0)) 



(-ir^i 



(2) 



with 6=1 — arccos( A) /tt, v is the spin-wave velocity and 
Ai,A2 are constants [7]. 

The non-oscillating term in Eq.(2) gives scaling dimen- 
sion for operator - d = 6/2 + 1/26 and from the com- 
mon formula [8] for mass gap m one has: 



1 



2-d 4-0-1/9 



(3) 



From Eq.(3) one could conclude that the magnetic field 
becomes irrelevant for A < — cos [tt-s/S] ~ —0.67 and 
the gap disappears for A < —0.67. This does not look 
physically reasonable, since the magnetic field destroys 
continuous symmetry of XX Z model and must produce 
the gap. In fact, due to nonzero conformal spin S* = 1 of 
the non-oscillating part of the operator it is necessary 
to consider higher-order effects in h [9]. The analysis 
shows, that in the pcrtiirbation series another critical 
exponent appears, giving for the mass gap 



7 = 



1-1 



(4) 



It turns out that the oscillating part of the operator 5^ 
gives another, more relevant index for the gap at A < 0. 
Let us reproduce usual conformal line of arguments for 
this oscillating part. 

The perturbed action for the model is 



S = So + h j dtdx S'^ix^t) 



(5) 



where 5o is the Gaussian action of XX Z model. Let us 
perform an infinitesimal renorm-group step with a scale 
factor A = 1 — so that x = Xx', t = Xt'. The cor- 
relation length changes as ^ = A^'. Then, the action 
becomes 



S' = So + hj d {Xt') d (Aa;') S''{Xx', Xt' 



Now let us estimate the large-distance contribution to 
the action of the oscillating part of the operator S^{x, t): 



h / dtdxS''{x,t) r^h dtY^ 



{-ly 



(„2 _ ^2^2)6/4 



hfdtY^ 



6n 



So, after rescaling we get 



h / dtdx- 



6x 



(a;2 -u2i2)|-M 



h j dtdxS''{x,t) /iA^-^/^y^ dt'dx' S''{x',t') 

and, therefore, the magnetic field scales as h' = hX^~^/'^. 
Expressing A as 



A=4=^=r^ 



f m \h 
we find that the mass gap is proportional to 

1 



T = 



i-e/2 



(6) 



Actually, the oscillating factor (—1)" in the correlator 
in some sense eliminates one singular integration over x, 
and into common conformal formula m ~ , where D 
is the dimension of space and d is the scaling dimension 
of perturbation operator, one should use D = 1 instead 
of conventional D = 2. 

The comparison of the expressions Eqs.(3), (4) and (6) 
shows that for < A < 1 the leading term is given by 
Eq.(3), while for -1 < A < by Eq.(6). Thus, one has: 



m ^ h 



< A < 1 
-1< A < 



(7) 



For example, rw^ h, when A ^ ±1 and m ~ h!^/^ for 
A = 0. 

In this respect the model (1) is different from the XX Z 
model in the staggered transverse field for which m 

^2/(4-9) < ^ [;^o]. 

The staggered magnetization (LRO) along Y axis be- 
haves as: 

{SD ~ {-l)^/e'^ - (-l)"m^/2 (8) 
Therefore, the LRO has also two different critical indices: 



\sy\) 



< A < 1 
1< A < 



(9) 



The critical indices v and r can be also found from the 
analysis of divergences of terms of perturbation series in 
h at N ^ oo {N is the system size). As it is shown [4] 
at A'' — > oo the ground state energy has a form 



(10) 



where v and r are given by Eqs.(3) and (6) and a, b are 
some constants. 

One can see from Eqs.(3,6) that — » 1 at A ^ 1 and 
r— >latA— >— 1. Hence, in both limits one of the sin- 
gular terms becomes proportional to h"^, and, therefore, 
gives a contribution to the susceptibility. It implies that 
in the symmetric points A = ±1 the susceptibility has a 
jump. For example, x = 1/4 at A = — 1 and x = 1/8 at 
A--l[4]. 
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The line A = 1. 

In the vicinity of the hne A = 1 it is convenient to 
rewrite the Hamiltonian (1) in the form 

H = Ho + V 

Ho = ^ S„ • S„+i — h ^ 
y = -5 ^5^5^+1 (11) 

where the parameter 5 = (1 — A) <C 1 is small. On 
the isotropic line A = 1 the model (1) is exactly solv- 
able by Bethe ansatz. The ground state of Hq remains 
spin- liquid one up to the transition point he — 2, where 
the phase transition of the Pokrovsky-Talapov type takes 
place and the ground state becomes completely ordered 
ferromagnetic state. Therefore, for /i < 2 and for small 
perturbation V we can use conformal estimations. 

The large distance asymptotic of the correlation func- 
tion on this line is 



(12) 



where a{h) is known function obtained from Bethe ansatz 
[11] and having the following limits: 



a{h) ~ 1 



21n(l//i) 



(13) 



and a(2) = 1/2. 

So, the scaling dimension of operator is rf,, = 
a{h)/2, and the scaling dimension of operator SfSfj^i 
is dzz = 4:dz = 2a{h). Since a{h) < 1, then the pertur- 
bation V is relevant and leads to the mass gap and the 
staggered magnetization along Y axis 



m ~ 5l/(2-rf„) ^ gl/(2-2a) 



{\sy\) 



(14) 
(15) 



Above consideration is valid also for the case A > 1 
(region 1) with the same exponents for the gap and the 
LRO. The only difference is that the staggered magneti- 
zation appears in this case along the Z axis. 

From the general expressions for the mass gap (3), in 
the limit ft, ^ 2 we obtain that m ^ g. 

The LRO in the vicinity of the point A = 1 , /i = 2 
vanishes on both lines: at A = 1 from (14) as 5^/*; and 

1 /8 

at h = he BS \hc — h\ ' (see Eq.(19)). Combining these 
facts we arrive at the following formula: 



(5^)^(-l)V/"|/^c-/i| 



1/8 



(16) 



which is in accordance with the exact expression for LRO 
on the classical line [3] . 

The behavior of the system near the point A = 1, h = 
is more complicated. As it follows from Eq.(7), for very 
small h the mass gap is m h, while on the other hand 
from Eq.(14) one obtains another scaling m ^ g 



ln{l/h) 



Really [12,4], there are two regions near this point with 
different behavior of the mass gap and a crossover line 
V5ln(l//i) ~ 1: 



m h. 



m ■ 



Jn{l/h) 



for V5ln(l//i) > 1 
for ^ln(l//i)«:l 



(17) 



The transition line h = ftc(A). 

Now we consider the behavior of the model in the vicin- 
ity of the transition line hc{A). For this we have used 
the Fermi-representation of (1). At first, in (1) we per- 
form a rotation of the spins around the Y axis by 7r/2 
(so that the magnetic field will be directed along the Z 
axis) followed by the Jordan- Wigner transformation to 
the Fermi operators a^,an- As a result, we obtain the 
Fermi Hamiltonian in the form 



HN N 



1 



1 + A 



2 

1 - ^ V- • , . + + 

H — 2^ sm k{a],aZk + a—kak) 



COS k)a^ak 



(18) 



Treating the Hamiltonian Hp in the mean-field approx- 
imation we find the ground state energy Eq and the ex- 
citation spectrum e{k). 

Main results following from the mean-field considera- 
tion are: 

1. The function e{k) has a minimum at fcmin, which is 
changed from 7r/2 at ft, = to zero at ft = ft-i(A) and 
kmi-a = for ft > fti(A). The gap in the spectrum £{k) 
vanishes at ftc(^) {he > h\) and for ft > fti is m ~ 
|ft — ftc|. The dependencies of fti( A) and ftc(A) are shown 
on Fig.l. There is the staggered magnetization along 
Y{Z) axis for A < 1(A > 1) at ft, < he and it behaves 
as ~ |ft — hc\^^^ at ft ^ he- The magnetization s = {S^) 
has a logarithmic singularity at ft — * ftc. 

These results show that the transition at ft = ftc(A) 
belongs to the universality class of the ITF model. 

2. The mean field approximation is exact on the classical 
line ft = ftci(A). 

3. In the vicinity of the point ft = 2, A = 1 the fermion 
density is small and the mean field approximation gives 
the accuracy in the energy, at least, up to or (2 — ft)*. 
For this case the gap is 



m 



[ft - ftc 



ft > fti 



m = 



2^/2 



he — h — 



32' 



h<hi 



(19) 



where ftc = 2— | — |2,ftl = ftc^fQ• 
It is interesting to compare Eqs.(19) with the confor- 
mal estimation of the gap m ^ g. The conformal ap- 
proach determines g dependence only, while Eq.(19) gives 
a prefactor depending on ft. 

The magnetic susceptibility x{h) is 
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X= — In 
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X 



ng \hc — h 

1 1 
\/27r \/hc - h' 



g > y/hc-h 
g < y/hc-h (20) 



As it follows froin (20) there is a crossover from square 
root to logarithmic divergence of x- 
The line A = — 1. 

On the line A = — 1 the model (1) reduces to the 

isotropic ferromagnet in a staggered magnetic field. This 
model is non-integrable, but it was suggested [13], that 
the system is governed by a c = 1 conformal field theory 
up to some critical value h = ho, where the phase tran- 
sition of the Kosterlitz-Thouless type takes place. In the 
vicinity of the line A = — 1 the Hamiltonian (1) becomes 

77 = -^S„.S„+i-/i^(-l)"5^ 

+ (1 + A) ^5^5^+1 (21) 

where (1 + A) <C 1 is a small parameter. 

It can be shown [4] that at hci{A) < /i <C 1 low energy 
states of the (21) is described by the XYZ Hamiltonian 



at ft, = ho, where P{ho) = 4. So, at the point h = ho 
the transition is of the Kosterlitz-Thouless type, and for 
h> ho the mass gap is exponentially small. 

Using the conformal invariancc of the model (21) at 
A = — 1 and h < ho we carried out finite-size calculation 
of the exponent /3(/i). The extrapolated function f3{h) 
agree well with the dependence n/h at ft <C 1 and /3 = 4 
at ho — 0.52. On the other hand, the mean-field approach 
gives rather crude value ho = hc{—l) = 0.69. 

In summary, we have studied the ID XXZ model in the 
transverse magnetic field. It is shown the spectrum of the 
model is gapped except critical lines on the (h, A) plane, 
where the LRO vanishes. Wc found the critical exponents 
of the gap and the LRO in the vicinity of these lines. 
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(22) 

The derivation of this mapping is based on the fact that 

the transition operator X](^l)"'5'n connects the low-lying 
states with the states with high energies ~ 2 only. The 
coincidence of the low-energy spectra of (21) and (22) 
for ft.ci(A) < ft, ^ 1 has been checked numerically. The 
spectrum of low-lying excitations of the XYZ model [14] 
as well as initial model (1) in the vicinity of the point 
A = — l,ft = can be described asymptotically exactly 
by the spin- wave theory [1] , which gives 



m = ftV(l + A)/2, A>-1 
m= v'(l + A)(l + A-/iV2), 



A< -1 



(23) 



The validity of the spin-wave approximation is quite 
natural because in the vicinity of the point A = — 1, ft, = 
the number of magnons forming the ground state is small. 

Wc note that the gap (23) for A > —1 coincides with 
the conformal theory result (7) and provides us with pre- 
exponential factor for the gap. 

On the line A = -1 the model (22) is the XXZ model 
and the correlation functions have the power-law decay. 
The scaling dimensions of operators Sf and Sf,Sf on 
this line are d,^ = /3(ft)/2 and dy = d,, = l/2/3(ft). The 
function /3(ft) is generally unknown, but at ft <C 1, where 
the mapping to XYZ model is valid, /3(ft) ^ ir/h. 

Strictly on the line A = — 1 at some vahie of ft, = ho 
the gap appears. It means that the magnetic field term 
is irrelevant at ft < fto (/?(^) > 4) and becomes marginal 
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